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ABSTRACT 

^ ■ We calculate the contribution to the ultraviolet background (UVB) from thermal 

\ emission from gas shock heated by cosmic structure formation. Our main calculation 

■ is based on an updated version of Press-Schechter theory. It is consistent with a more 

empirical estimate based on the observed properties of galaxies and the observed cos- 

, mic star formation history. Thermal UVB emission is characterized by a hard spectrum 

extending well beyond 4 Ry. The bulk of the radiation is produced by objects in the 
mass range 1O^^~^^M0, i.e. large galaxies and small groups. We compute a composite 
UVB spectrum due to QSO, stellar and thermal components. The ratio of the UVB 
• intensities at the H and He Lyman limits increases from 60 at z = 2 to more than 300 

P-i' at z = 6. A comparison of the resulting photoionization rates to the observed Gunn- 

Q . Peterson effect at high redshifts constrains the escape fraction of ionizing photons 

^ ' from galaxies to be less than a few percent. Near 1 Ry, thermal and stellar emission 

I are comparable amounting to about 10 %, 20 % and 35 % of the total flux at redshifts 

■ of 3, 4.5 and higher, respectively. However, near the ionization threshold for He ii , the 

thermal contribution is much stronger. It is comparable to the QSO intensity already 
at rcdshift ^ 3 and dominates at redshifts above 4. Thermal photons alone are enough 
to produce and sustain He ii reionization already at z w 6. We discuss the possible 

^ I implications of our results for the thermal history of the intergalactic medium, in 

particular for He ii reionization. 

Key words: cosmology: large-scale structure of universe — radiation mechanism: 
thermal — shock waves 



1 INTRODUCTION 

Neutral hydrogen in the intergalactic medium (IGM) produces a forest of resonant Lya absorption lines in the spectra of 
high-redshift quasars. The connection of these featur es to the structure formation process has now been now firmly e stablished 
using N-body/hydrodynamic numerical simulations llCen et al.lll994t[zhang et al.lll995l : iMiralda-Escude et al.lll993) . There is 
consensus that the observed IGM temperature results from a balance between photo-ionization heating and adiabatic cooling 
due to the Hubble expansion. Such photo-heating is provided by the extragalactic ultraviolet background (UVB) whose nature, 
origin and evolution have, therefore, been subject to considerable investigation. 

At low redshifts, the ionization balance is consistent with a pure power-law ionizing spectrum. Traditionally QSOs have 
considered as the mai n sources of ionizi ng photon. However, a number of recent theoretical and observational developments 
are at odds with this. iKim et all ll200lll find that the break in the redshift evolution of absorbers occurs at lower redshifts 
than predicted by numerical simulations usi ng a standard QSO ionizing background (Haard t fc Madau 199S,) , hinting at an 
incomplete description of the UVB. This led lBianchi et all i200ll . BOl hereafter) to recompute the UVB as a superposition of 
contributions from QSOs an d galaxies. Those authors adopted an escap e fraction of ionizing photons from ga laxies, /esc — 10%. 
In spite of many theoretical JPove et alJbOOd: IWood fc Loebl20 00: H aehnelt et al.ll200ll:ICiardi et alJl2002^ and observational 



* fm@MPA-Garching.MPG.DE 



2 Miniati et al. 



studies jLeithe rer et alj |l995l : iHeckman et aJboOlt ISteidel et alJlioOlh the determination of this param eter remains hiRhly 
uncertain. In i lT^ bclo w we use recent measurement s of the Gunn- Peterson effect in high redshift quasars jBecker et aljliooil : 
[McDonald & Mir alda-Escudell200ll : iFan et ;il]l2002l ) to argue for values of /esc as small as a few % for high redshift galaxies 
( tl3.3L in agreement with independent estimates by ^Fcrnandez-Soto ct al. (2003). These authors set a 3(t (statistical) upper 
limit /esc ~ 4% using photometry of 27 spectroscopically identified galaxies with 1.9 < z < 3.5 in the Hubble Deep Field. 

Sim ilarly, high resolution si mulations using a standard cold dark matter model and a standard QSOs ionizing back- 
ground iHaarc|t fc Madad Il99fil). , produ ce a Ly-a forest lines with a minimum width significantly below that observed 
jTheuns et al. "l998a': 'Brvan ct alJ ll999^ . This appear s to re quire additional heat sources, for example photo-electric dust 
heating I ^Nath ct al. 1999), r adiativ e transfer effects iAbel fc Haehnelt 19991) . or Compton heating by X-ray background 
photons |"Madau fc Efstathioui ri99fll) . 

The situation around redshift z ~ 3 is more complicated. Apparently the spectrum must be very soft, with a large break 
at the He"^ edge. Songaila ( 1998) reports an abrupt change of the C iv /Si iv ratio at z « 3. This may indicate a rapid 
and significant change in the shape of the ionizing spectrum, perhaps due t o He ii reion i zation . This is corroborated by the 
detection of patchy He ii Ly— a absorption at similar redshifts. In addition, iHeap et al.l i200(i) analyzed the Gunn- Peterson 
He II absorption trough at z — 3.05 found the UVB to be characterized by a high soft ness parameter^ S « 800, in contrast 
to the harder spectrum deduced at z = 2.87. However, a more recent VLT/UVES study ("Kim et aljr2002ri using 7 QSOs finds 
no strong discontinuity for the quantity C iv /Si iv around z — 3 and suggests that it might not be a good indicator of 
He II reionization. These results are in agreement with the analysis of feoksenberg et alj (12003) who studied the properties of 
met al absorption system s between redshift 1.6 and 4.4 in the spectra of 9 Keck/HIRES QSOs. 

ISchave et alJ (l200(i) concluded that the IGM t empe rature evolution differs considerably from simple expectations based on 
QSOs as sole ionizing power input. iRicotti et al.l (|200(J) found a similar result with an analogous but independent approach. 
Both analyses suggest at z ~ 3 the IGM temperature (of the gas at mean density) undergoes a "sudden" jump which is 
interpreted as associated with He ii reionization. It is worth pointing out that lRicott i et alj ||2000|) re cover higher temperature 
values and a quite smoother temperature jump than lSchave et al] i2000l) . whereaT lMcDonaMetS^ (1200 ll) using similar data 
to Shaye et al. and a more conservative analysis did not find the same temperature step at redshift 3 (for IGM gas overdense 
by a factor 1.4 with respect to the mean). This provides a measure of the possible uncertainties in these experiments. 

In this paper we explore another source of UVB ionizing photons, namely thermal emission from shock-heated gas in 
collapsed cosmic structures. We find that thermal radiation, produced mainly in halos with temperatures between 10^ K and 
a few X 10^ K, is characterized by a hard spectrum with many photons above the H i and He ii ionization thresholds. As noted 
above we limit the escape fraction of UV ionizing photons from high redshift galaxies, /esc, to a few %. Thermal emission 
provides a significant fraction of H i ionizing photons at redshift J5 3. In addition, we find that thermal emission plays a major 
role in the reionization of He ii , being comparable to the QSOs at redshift ~ 3 and dominating the flux at z > 4. In fact it 
turns out that thermal photons alone are enough to cause He ii reionization at 2: ~ 6. Our study is based on well understood 
emission mechanisms, bremsstrahlung and line emission from optically thin thermal plasma, so we expect our results to be 
robust. 

Thermal emission is not the only possible source of ionizing photons in collapsing structures. An alte rnative is inverse 
Compton emission by relativistic electron s accelerated i n large scale structure shocks jMiniati et ahlliooJ) . However, based 
on methods similar to those pres ented in | Min iati| l|2pp2j) we find such a component to be negligible compared to that from 
QSOs, a result in agreement with iRandall fc SarazinI J200J) . 

The paper is organized as follows. In f|5|we present the details of our model, that is we describe how we solve the radiative 
transfer equation and estimate the radiation due to QSOs, stars and, particularly, thermal emission. Our results are presented 
in S|H]where the relevant features of the transmitted radiation fiux are described. Finally, in S0]we discuss the thermal evolution 
of the IGM and in 33 we summarize our main results. 



2 MODEL 

2.1 Radiative Transfer 

The mean specific ionizing flux, J(i/o,Z o), observed at frequency, Vo, and redshift, Zo, is the solution to the cosmological 
radiative transfer equation which reads l|Peeblejll993ll 

T( ^ " r -r.,,(.„,.„,z) [n^(l + ^)^+nA] "'^' ( l + Z, \^ 

where c is the speed of light, Hq is the Hubble parameter, Q,m and J^a are the matter and vacuum energy densities respectively^, 
and j{y,z) is the volume- averaged proper spectral emissivity computed at emission redshift, z, and at the appropriately 
blueshifted photon frequency u = Vo + 2)/(l + Zo). In addition, Teff{uo, Zo, z) is the effective optical depth at frequency Uo 



^ Note that S is defined as the ratio of the photoionization rate in H I over the one in He II ( Ai.'Il . 

^ We assume a flat cosmological model with normalized Hubble constant hyo = Ho/70 km s^^ Mpc~^ = 1, Qm = 0.3, Ha = 0.7, 
as = 0.9. 
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due to absorption of residual neutral gas in the IGM between Zo and z. Following IParesce et al.l il98Ct) . for a distribution of 
discrete absorbers in an otherwise transparent (i.e. ionized) medium we write 

poo 

r,a{uo,Zo,z)= dz' dNm /{N-^^ ^ , z'){l ~ e'^^-"'^). (2) 

J Zo Jo 

The distribution of absorbers as a function of the H i column density and redshift, f{Nm,z') = N/dNmdz' , is typically 
derived from counts of Lya lines in QSOs absorption spectra. It is assumed that the number density of absorbers evolves with 
redsh ift as dN/dz' cx (1 + implying an evolution of the effective optical depth of the Lya forest as r^*^ ~ (1 + z)~' iZud 

Il993^. This was the approach taken in BOl who, based on observations of the Lya forest in the redshift range 1.5 < z < 4 
jKim et aLlbOOlTl . determined a value of 7 ~ 3.4. For redshifts higher than 4, however, recent observations, including the 
detection of a nearly complete ab sorption through in the QSOs spe ctra discovered by the SDSS at z ~ 6, imply a much 
stronger evolution of that function iBecker et al ]|200ll : iFan et alJl2002ll . Therefore, for z > 4 we retain the power-law behavior 
fo r t'^^ but d eterm ine the power-law index 7 through a comparison of the power-law models with the new measurements 
of iFan et al.' ('2002"). In Fig.H] we report the power-l aw mode ls for different values of 7 together with observed data points 
taken from McDonald & Miralda-Escudc (200j) and lFan et a l. (2002). Actually the data do not exhibit a smooth behavior 
and, therefore, cannot be described in full detail by any of the smooth power-laws. The features in the data are likely due to 
inhomogeneities in the r eionization process or to a drop in the ionizing flux in the aftermath of cosmological reionization (e.g. 
ICen fc McDonaldll2002l) . In addition, as we approach 2; ~ 6, the data indicate a sudden increase of the IGM optical depth. 
Both of these features in the redshift dependence of are reflected in the redshift evolution of the ionization rates that 
[McDonald & Miralda-Escudc (2001) and Fan et al. (2002) inferred by modeling the distribution of absorbers in the IGM in a 
fashion similar to ours. Here we adopt a single power-law with 7 — 5.5 at z > 4, which a pproximately describes the behavior 
of T^^iz). Howe ver, it is part of our objectives to compare the measured ionization rates iMcDonald fc Miralda-Escudell200ll : 
iFan et al]l2002ll with those produced by the radiation processes explored in this paper. We will heretofore need to bear these 
approximations in mind in order to properly interpret the results. 

Each absorber is characterized by an optical depth r(z/') = A^hio"hi(j^') + AfHciio"Hcii(j^')i due to H i and He 11 ionization 

I He I ionization is neglig ible). The value of A'^hoII can be derived from A''hi by studying the radiative transfer within each cloud 
Haardt fc Madaulll996f) . For simplicity, we have assu med that all c louds are optically thin at the He 11 ionization threshold, 
which yields A^Heii/A^Hi ~ 1.8 J(13.6eV)/ J(54.4eV) (lMadaulll99ih . While this approximation slightly underestimates the 
optical depth at hu > 54. 4 eV, it does not sever ely affect the opacity of the most abundant (A^hi < lO^'^ cm"'^) clouds for 
hu < 100 eV (cf. Fig. 2 of lHaardt fc Madaulll99d) . 

In order to compute the mean specific ionizing fiux eq. (0 is solved iteratively, as required by the dependence of Nhcu/Nhi 
on z). For the total emissivity we take j = jq + jc + Jt, that is the sum of QSO, stellar and thermal components. In the 
following sections we describe how each component is computed. 



2.2 Radiation from QSOs and stars in galaxies 

The contributions of QSOs and stars to the ionizing UV background adopted in this paper are similar to those derived in 
BOl. We briefly summarize here the m ain assumpti o ns of the model. The QSO emissivity assumes a lumi nosity function tha t 
follows the double power-law model of iBovle et al.l lll988ll . For 2 ^ 3, we adopt the parameters given in lBovle et alJ ll200(J) . 
ob tained by fitting a sample of over 6000 QSOs with 0.35 < z < 2.3. At z > 3, we include the exponential decline suggested 
bv lFan et alJ (1200111 . describing the dramatic reduction in QSOs number density at high red shift. Finally the QSO spectrum 



in the ionizing UV range is modeled as a simple power law, oc ly' (IZheng et al.lll997fl 



For the stellar component, we assume a star formation rate that is constant from high redshifts to z ~ 1, then rapidly 
decreases to local values, as indicated by several galaxy surveys in the rest-frame no n-ionizing UV (iMadau. et al.l Il998l : 
ISteidel et alj [19991) . Synthetic galactic spectra (produced with the 2001 version of the iBruzual A. fc Chariot! Il993l code) 
have then been used to calculate the emissivity for the ionizing UV as a function of z. Finally, the internal absorption of 
radiation by the galaxy interstellar medium was modeled by a redshift-independent value for /esc, the fraction of Ly-continuum 
photons that can escape into t he I GM. As w e shall show in [0 the values for the ionization rates recently measured by 
McD onald & Miralda-Escude ( 200 J) and lFan"e t al. (2002) constrain /osc to a level of a few %. In the following we will usually 
set /esc ~ 1%. Wc sho uld also mention fo r completeness that we have neglected contributions from radiative recombinations 
in the absorbers ( Haardt fc Madciulll996l) . It would be possible to correct for this omission, but it is not strictly necessary for 
comparing the contributions from the various emission processes explored in this paper. 



2.3 Thermal Emission 

In this section we compute the contribution to the ionizing background flux from thermal emission from gas accreting onto 
dark matter halos. Since we calculate this for the flrst time here, we outline the assumptions of the model in more detail than 
for the QSOs and stellar components ( i\'2.'2t . 

The gravitational c ollapse of matter d ensity enhancements during structure formation drives supersonic flows which 
eventually shock (e.g. iMiniati et al.l liood) . The shocked gas is coUisionally ionized and heated to temperatures T ~ 
10*^ K («/100 km s"^)2 so that dense regions can efficiently radiate away thermal energy through bremsstrahlung and line 
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Figure 1. Multi power-law fits of the redshift evolution of r""' oc (1 + z)'' . Up to 2 ^ 4 a single curve corresponding to 7 ~ 3.4 fits 
the data. For higher redshifts we sh ow t he cases of 7 = 5 (dot), 5.5 (dash) and 6 (long dash) respectively. Data point are taken from 
iMcDonald k. Miralda-Escudj ^2nnl^ and lFan etlj] yfl 



emission JWhite &: Reej|l978h . At these temperatures and where the gas overdensity is above a few hundred collisions ensure 
the fractions, nHi/nmi ~ 10"'^ and riHcii/nHciii ~ 10~^, for neutral hydrogen and singly ionized helium respectively. It can be 
shown that under these circumstances the column density associated with the gas within the virial radius produces negligible 
absorption both at 1 and 4 Ry (see Appendix^. However, for cooler structures collisional ionization is less efflcient and 
the gas becomes optically thick for virialization temperatures T ~ 10^ K (see Appendix 0. Most of the thermal emission, 
however, is produced in objects with T > 10® K. . Therefore, in the following we will assume that thermal radiation is 
characterized by an escape fraction of order 1. 

In the following, the description of the evolution of the baryonic gas, responsible for the emission of th ermal radiation, 
is rather simplified . Althoug h in principle it could be pursued in some detail through semi-analytic schemes jWhite fc FrenkI 
ll99ll:lKauffmann et al. 1993: Ivan den Boschl2'ooll . the latter provide considerably more information than we require and their 
level of sophistication is well beyond the scope of our current study. However, we do test our model predictions for the amount 
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of cooling that must have occurred against independent empirical estimates. These are based on the observed cosmic star 
formation history and the distribution of stellar mass as a function of halo masses, that we infer from recent SDSS data. The 
comparison is detailed in tl2.4l In addition, we comment on possible additional effects neglected in our simplified approach. In 
general, however, the agreement that we find between our predictions and our empirical estimates is encouraging and supports 
our adopted approach. 



2.3.1 Collapsed Halos 

The volume averaged comoving thermal emissivity in units of erg s~^ 



cm ^ Hz ^ can be written as 



iT(!^,z) 



dM —{M,z) Me[T{M),u] min 



TcooliM, z) 
THubbleiz) 



(3) 



where dn/dM is the comoving number density of collapsed dark matter halos of mass M at redshift, z, and e{T,i^) is the 
spectral emissivity per unit mass. The last term is the ratio of the cooling time to the Hubble time at redshift z and accounts 
for the fact that any parcel of gas em its only as long as i s allow ed by its supply of therm al energy. 

To describe dn/dM we adopt the lPress fc Schechterl (Il974ll formalism as updated bv lSheth fc TormenI lll999f) . This gives 
the mass function 



dn 
dM 



iM,z)^A{-) 



Po 

M2 



-2p\ 



dlna{M) 



dlnM 



exp 



(4) 



where po is the current (2 = 0) average mass density of the universe, £> = a^^^Sc{z)/o-{M), a[M) is the mass variance of linear 
density perturbation corresponding to mass M at the current epoch and 5c (2) is the critical linear overdensity e valuated at 
present for a spherical perturbation that collapses at reds hift 2. Both a(M) and Sd z) are taken from lKitavama fc~Suto ( 199'^ 
and the parameters {A, a,p) = (0.322, 0.707, 0.3) are as in lSheth fc TormenI lll999f) . According to the spherical collapse model, 
for 2 S> 1 the virialized dark matter halos are characterized by an overdensity Ac — IStt^ with respect to the critical value, 
per. The baryonic gas is assu med to settle in approximate hydrostatic equilibrium within each halo at its virial temperature 
(e.g. iBrvan fc Normanlll998tl 



T(M) = 1.4 X lO^K h. 



2/3 



/ M 



2/3 



V0.3J Ustt^J V 5 ) 



(5) 



and with a density pi, = ft, A^pcr, where fi, is the halo baryonic fraction. 

The spectral thermal emissivity, e(T,iy), is computed through the code bv iRavmond fc SmithI dl2Z2l version 1992)^ 
where as the cooling function, A(r, Z), necessary to compute the cooling time, was taken from the work of lSutherland fc Dopital 
il993h . Our choices are appropriate because, for regions characterized by an overdensity of order 100 or greater at redshift 
~ 2 — 3 and higher, the ionization equilibrium is primarily determined by plasma colli sions. We adopt an ave rage metallicity 
within the collapsed halos Z = Zq/20, thought to be typical for high redshift objects iMartinelli et al.ll200Cf) . 

Contributions to the integral in eq. Q are limited below a certain mass threshold. For halo masses below M;c = 
2.6 X 10" (!r2m/0.3)-^/^ /(I + 2)]^/^ M0, where £i3.6ev is the photon energy in units of 13.6 eV, emission of 

thermal radiation is depressed exponentially because the low virial temperature is insufficient to ionize hydrogen. The most 
stringent constraint for low mass objects is due, however, to cooling. In fact, we find that the cooling correction in eq. ^ is 



'^cool 
THubble 



^^^^JjjllAble_ 

ngasA{T,Z) 



: a0.4h: 



-1/3 



/ M 



\ 1013 Mq 



2/3 



-1/2 



VISttS 



-2/3 



A(T, Z = Ze/20) 



10- 



^erg cm-^s" 



' / fb 
\0.3j \0.15j ' 



(6) 



where Ugas is the gas number density, Z is the gas metallicity, A is the cooling function and fcs is Boltzmann's constant. For Z ~ 
ZfTi/2 0, the normalization for A given in eq. ||SJ holds to a good approximation for 10^'^ ^ T 5; 10^ (cf. [Sutherland fc Dooit^ 
Il993fl . If the gas contraction induced by radiative cooling occurs isobarically the amount of radiated energy equals the enthalpy 
of the system and, a = ygaa/i'Jgas — 1) = 5/2. We point out that the halo life-time, instead of the Hubble time, should perhaps 
be used in eq. Q, but we will neglect these details in the current investigation. Thus, the above scenario should provide at 
least a lower limit to the potential amount of energy to be radiated (even without considering the effects of feedback, see 
below) . 



2.3.2 Cooling & Feedback 

It is well known that feedback of energy and momentum associated with formation of stars plays a fundamental role in 
regulating the dynamics of the IGM, although a clear and coherent understanding of how is still lacking. For the purposes 

3 The code assumes the emitting plasma in collisional equilibrium and neglects the effects of an ambient radiation field. It includes the 
following atomic processes: collisional ionization, collisional excitation followed by auto-ionization, radiative recombination, and dielec- 
tronic reco mbination. Collisional ioni zation and excitation are assumed to be produced by free electrons with a Maxwellian distribution 
of energies <Ravmond fc Smithlll977t) . 
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of the present investigation, the effects of feedback can be taken into account by computing the increase in the coohng time, 
Tcooi, caused by the injection of additional energy. The average energy deposited in each halo of mass M can be computed 
through two quantities: the fraction /, of the baryons (fbM) that is converted into stars and, the amount of energy, (^sn, 
released by supernova explosions per unit of stellar mass formed. When normalized to the halo volume-integrated enthalpy, 
W = {5/2)nkBTV — {5/2) fbM ksT/ unip, where V is the halo volume and fimp is the gas mean molecular weight, the 
deposited energy is 

A^FB ^ f.fbMCsN ^ (J^\ ( CSN \ 

W W ■ V0.15y' I^1049erg MqV 

The value assumed for (^sn varies in the literature. If we take a ratio of supernova per total mass converted into stars 
~ 4 X 10""^, as sugge sted by a s tandard Salpeter initial mass function, then Csn ^ 4 x lO'^^ erg M^^, although base d on 
metallicity arguments ISilkl ((So^) suggests a value 4 times larger (see also lBookbinder et alJlToSol : ICen fc OstrikeI^ll992^ . We 
notice, incidentally, that the deposition of the assu med amount of ener gy corresponds to fJ^mpft,(sN ~ 1 keV/part. roughly 
as observed in the core of small groups of galaxies dPonman et al.|[l999^ ■ Obviously, for the assumed parameters, the energy 
injection will only substantially affect gas within halos with temperatures ^ f ew x 10^ K. This result is in agreement with the 
more sophisticated feedback prescriptions presented in IWhite fc FreiiS dlPgiTl . 

With the additional energy expressed in eq. ^ the radiative lifetime of a halo, as in eq. Q, will be extended by a factor 
(1 + tjAEfb/W). Here we have introduced yet another parameter, r;(M), the fraction of energy deposited through feedback 
processes that is available for conversion into thermal radiation. We do so in order to account for various effects produced 
by feed back. For example, when AEfb/W ^ 1 gas is more likely to be blown out of its host halo, thus strongly inhibiting 
cooling jMacLow fc FerrarallTgggI) . This would imply rj ^ 0. In addition, the e volution of a cooling core will be altered and 
possibly disrupted by the occurrence of merger events (cf. IWhite fc Reeslll978f) . making 77 < 1. However, as we shall show in 
the following, most of the radiation is contributed by halos within a narrow mass rang e. Thus, in practi ce, we only explore 
the case 77 ~ 1/2, which we consider appropriate for the mass range of interest here (cf. iMori et alJl2002|) . 



T{M) 



(7) 



2.3.3 Halo Emissivity Distribution 

In Fig. 121 we show which halos are primarily responsible for the production of thermal far UV emission. For this purpose 
we define arbitrarily a UV luminosity LuviT), as the integral between 13.6 and 100 eV of jT(v,z) given in eq. We then 
plot for different redshifts and as a function of the temperature of the emitting halo population, histograms of dL„„/dlogr, 
that is the total UV luminosity of halos with virial temperature in a logarithmic interval centered on T. According to our 
model, most of the emission is produced by halos with temperatures between 10® K and a few xlO^ K, corresponding to 
masses 10^^~^^ Mq. The contribution from halos below this range is reduced because of their short cooling time. This may 
be seen by comparing with the highest (no-cooling) curve which corresponds to the case in which the corrective term that 
accounts for the effects of cooling is omitted. Notice that the "bumpy" shape of these curves for r 10^ K is a reflection of 
the structured temperature dependence of the cooling function. On the other hand, the bend at the high mass end is solely 
due to the paucity of massive halos above a certain threshold which depends on both redshift and a{M): for a{M) oc , 
dLuv/dlogT oc dLuv/dlog M oc M'^^^'^^^. The line labeled 'w-cooling-|-feedback' corresponds to the feedback case described 
in Sec. 12.3.^ As expected, the extra injected energy primarily affects small structures which are now able to produce a larger 
amount of radiation (by reradiating the injected energy). Compared to the no- feedback case (no- feedback) the spectra in Fig. 
|l|are increased by a factor of a few in the low energy UV range but are basically unchanged in soft X-rays. This is a reflection 
of the fact that feedback mostly affects smaller (colder) halos. In any case the objects which dominate the thermal emission 
are only weakly affected by the specific feedback prescription we adopt. 



2.4 Comparison with Observational Data 

We will now attempt to test our model against observ able quantities that are related to the occurrence of radiative cooling. 
For this purpose, we have used SDSS data presented in lKauffmann et al.l (|20o3) for a representative sample of more than 10'' 
galaxies, to construct a function, g*(M), describing how the total stellar mass is distributed over halo mass in the present 
day universe. The dataset includes the following quantities of interest to us: colors in the g, r, i bands, stellar masses and 
concentration index. 

For each galaxy, the halo mass is inferred from its luminosity through either the TuUy-Fisher or Faber- Jackson relations 
depending on whether the galaxy is disk or bulge dominated, respectively. We first divide the sample into spirals and ellipticals 
on the basis of the concentration index defined as the ratio of the Petrosian half-light radius to the Petrosian 90% light radius, 
C = rpo/rsQ. Thus, a galaxy is classified as spiral or elliptical/spheroidal depending on whether or not C > 3 ( Shimasaki^^^ 
|20^j). The SDSS color indexes are converted to Johnson's UBVRcIc system according to the relations given bv Smith et al. 
i2002l) . Next, for spirals we use Verheijen's (2001) relation between the /—band magnitude. Mi, and the rotational velocity 
Vfiat, measured in the outer parts of the galaxy disk where the rotation curve flat tens out. For ellipti cals we employ the 
relation between B— band magnitude, Mb, and dispersion velocity as determined bv lBender et al.| lll99d) . Temperature and 
velocity are then related according to ksT = (l/2)/impii^ and temperature and mass as in eq. The results for (?, are 
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Figure 2. Histograms of the integrated thermal UV emission between 13.6 and 100 eV as a function of the virial temperature of the 
emitting halo population for four different redshifts. The various lines correspond to cases with (w-cooling) and without (no-cooling) 
accounting for the effects of radiative cooling, and with cooling+feedback ( w-cooling+feedback) . The thick line (SDSS) was obtained 
based on the observed cosmic star formation history and the distribution of stellar mass as a function of halo virial temperature as 
reconstructed from SDSS data (see Sec. l2.4l for details). 



shown in Fig.|21 where the histogram represents the fraction of stellar mass in the local universe as a function of the host halo 
virial temperature and mass. 

The formation of a mass Af« of stars implies that the thermal energy /enthalpy associated with the equivalent gas mass, 
aM^kT, must have been radiated away. Here T{M) is the virial temperature of the host halo and a — 5/2 as in eq. ®. In 
fact, this is just a lower limit. More precisely, the amount of radiated thermal energy is connected with the rate at which 
cold gas forms inside a halo, Mcoid, rather than the rate at which stars form, M*. Obviously, averaged over the halo lifetime 
Mcoid ^ Ma,; however, there is no straightforward way of inferring the former from the latter. A number of empirical facts 
and theoretical arguments suggest a proportionality between the cosmic-volume averages of the two quantities. 
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Figure 3. Histograms representing g«, i.e. the fraction of stars in the local universe hosted in halos of virial temper ature T, and mass 
A/. M and T are related through eq. where we set z = 0. The histogram was built on the basis of the SDSS data jKauffmann et all 



{Mcoid) = A(Af, z){KU){M, z) (8) 

with the proportionality factor A(M, taking values of order of several at high redshifts (in Appendix^ we provide an 
estimate of A(M, z) based on Press-Schechter formalism). Since we only have poor statistical information about how the time 
history of star formation depends on the halo mass, we assume that, on average, stars in halos of a given mass M formed at 
the measured cosmic star formation rate (see, e.g. BOl) multiplied by our stellar mass distribution function, g, (T). That is 

( Af, > (M, z) = <?, (M) SFR(z) . (9) 

This is the simplest assumption compatible with the results in Fig. |3 Thus, the differential amount of energy radiated as a 
consequence of star formation by halos of mass M is 
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which amounts to a differential volume-averaged thermal emissivity 

The last term in the above equation describes the spectral distribution of the radiated energy. Using the above 'empiri- 
cal' thermal emissivity, and the cosmic star formation rate SFR{z) summarized in BOl, we have recomputed the quantity 
dLuv/dlogT presented in Fig. |5]in order to compare it with our model predictions. The results are illustrated by the thick 
(SDSS) line plotted in Fig.|5|for each redshift (the adopted value of A was computed as detailed in Appendix f5t. 

The comparison between model and 'empirical' predictions is only meaningful for halos with TcooI < THubUe (no star 
formation can occur otherwise), i.e. with virial temperatures somewhat below the overlap point of the curves labeled as 
'with-cooling' (cooling accounted for) and 'no-cooling' (no cooling) in Fig. |21 The actual dividing line should occur a few 
times below this point of overlap, typically above several xlO^ K. This is because the cooling suppression factor in the 
former curve is based on the assumption that the halo age equals the Hubble time and is thus a bit overestimated. With 
this clarification, we conclude that the 'empirical' curve agrees quite well with the model predictions for star forming halos 
with virial temperatures above a few xlO^ K, where the cooling effects are strong. Below xlO^ K the agreement worsens 
considerably. This is somewhat expected and is tipically attributed to the suppression of star formation and/or gas blow-away 
due to feedback effects. Nevertheless, this is not too worrisome since the contribution of these objects to the UV background 
is not crucial, even in the 'feedback-case'. 

Finally, we note that the model and empirical curves have been derived using completely different methods and assump- 
tions, with values for the few free parameters involved (e.g. 6) taken from the studies in which these were introduced, rather 
than in order to improve our results. Nevertheless, the two approaches give very similar estimates of the UV thermal emission; 
this is encouraging as, in principle, the agreement could have been much worse. 



3 RESULTS 

3.1 Ionizing Flux 

Fig. 0]shows the mean ionizing flux due to thermal emission with (?7 = 0.5) and without feedback effects, that due to QSOs, 
and that due to galaxies (/esc = 1%) at four different redshifts. The imprint of the H i and He ii IGM absorption, with a 
significant reduction of the flux close to the Lyman limits of the two species, is clearly visible in each spectrum. We flnd that 
thermal emission provides an important contribution to the average UV radiation flux. Apart from this, the three sources 
of radiation have very different spectra and redshift dependence. At z ~ 3 and near the H i ionization threshold, thermal 
emission with feedback effects corresponds to ^ 10 % of the QSO contribution and is comparable to the stellar component. 
When feedback effects are not included the thermal flux in this spectral region is reduced by a factor of a few. At higher 
redshifts the thermal flux becomes progressively more important relative to the QSO component, whereas its relation to the 
stellar component is virtually unchanged. This latter feature makes sense because the production of stellar photons is linked 
to the formation of stars and, therefore, of cold gas, which in turn produces the thermal photons. Since the two components 
were modeled independently, this provides a nice consistency check. At photon energies near the He ii ionization threshold 
the trends are similar, except that now the thermal flux plays a more dominant role because of its relatively hard spectrum. 
Thus, in this energy range, and independently of feedback effects, thermal emission dominates the stellar contribution (even 
for /esc ^ 10%). It is comparable to the QSOs flux at 2 3, and becomes the dominant source at redshifts ^ 4. 



3.2 Softness parameter 

Both for the individual spectral components and for the total spectrum, we have computed the redshift evolution of the 
spectral softness parameter, 

Sl^^, (12) 

JHcII 

where Jhi and JhcII are the UVB intensities at the H i and He 11 Lyman limits respectively. These are shown in Fig. |^ As 
already anticipated, the thermal emission is characterized by small values of Sl, ranging from a few at low redshift to 100, 
at 2 « 6, with minor differences introduced by the adopted feedback prescription. A similar evolution is seen for the QSO 
spectra, although with higher Sl values in the range 80-300; the stellar component instead shows less pronounced evolution, 
increasing by only a factor of ~ 3, but maximal Sl values that exceed 1000 at 2: 3. The composite spectra (feedback case) 
Sl evolution resembles very closely that of QSOs, even at high redshift: this is somewhat fortuitous as at z 4 the UVB 
is dominated by t he sum of stellar and thermal contribution. The above values can be compared with the available data. 
iHeap et al.l (I2OOOI) . from an analysis of the quasar Q0302-003 (2 = 3.286), flnd that around 2 = 3.2 the ratio of the H i to 
He II photoionization rates, is larger than 400. As Sl ~ FHi/rHcii, where r depends on the spectrum shape and is about 
4 for our composite spectrum, from Fig. |^ we find that this value corresponds to our estimate for 2 = 3.4, a very close match. 
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Figure 4. Mean ionizing flux due to: thermal emission with (higher, soHd) and without (lower, solid) feedback effects, QSOs (dash) and 
stellar (assuming f^sc = 1%, long-dash) for four different redshifts. From left to right, the vertical arrows close to the horizontal axis 
indicate ionization thresholds of H I , He I and He II respectively. 



In our model the results of Heap et al. would not be necessarily caused by a jump in the parameter S at z ~ 3, but by the 
evolution of the composite spectrum, particularly the decrease in the ratio of thermal to the QSO flux. In the final Section 
we will reconsider this point. 



3.3 Photo-ionization Rates 

Fig. Elshows the photo-ionization rates in units of 10~^^ s~^ defined as 

ri2(.) ^ TT^ffirr = 4. r -4-) 4^ (13) 
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Figure 5. Redshift evolution of the softness ratio Sl for different emission cases: thermal, no feedback (dot); thermal, with feedback 
(thin solid); stellar (long dash); QSOs (short dash); total, no feedback (higher thick solid); total, with feedback (lower thick solid). 



for the various UV radiation components discussed above and as a function of redshift with the same notation as in Fig. 01 
In a ddition, wo also plot the photo-ionization rates inferred from the observed Gunn-Peterson effect in high-z QSO spectra 
^McDonald & Miralda-Escudc 2001; Fan et al. 2002). 

The left panel of Fig.|S]is relative to ionization of H i and contains a number of important features. First, we notice that 
although the emission from QSOs is able to produce the ionizing flux observed at 2 ~ 2 — 3, it falls short at higher redshifts, a 
well known fact. In o ur formulation, it results from the assumed ra pid decline of the QSO nu mber density for 2: > 3, as derived 
from the SDSS ( H2.2lFan et al.ll200j) . We note here, however, that iMeiksin fc Whitd i2003l) in a recent paper questioned this 
view: by tweaking the bright end of the SDSS QSO luminosity function and allowing for objects fainter than the survey limit, 
QSOs alone can in principle account for the UV background at z > 4. 

As a second result, the comparison between the meadured values of ri9 lMcDonald fc Miralda-Escudl i200j) : lFan et all 

i2002h and the stellar ionization rates assuming fesc = 1% (dash curve) imply, according to our model, that fesc has to be 
smaller than a few %. In fact, had we assumed, e.g. , f^sc — 10% the observed F would have been overpredicted by a factor 5 
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Figure 6. Left: Photo-ionization rates defined in ea. ll3l as a function of redshift for ionization of H I due to emission from: QSO (dash), 
stellar for fesc = 1% (long dash) and shocked IGM with feedback effects (solid) and without (dot). The thi ck solid line is the total consid- 
ering the f eedback ca s e (or, alternatively, the no-feedback case and feac = 2%). The data points are from lMcDonald fc Miralda-Escudj 
i20oil) and lFan et'al] 12002^ after correction for our cosmological model. Right: same as for left panel but now for coefficients relative to 
He II (right). The stellar component (long dash) has now been multiplied by a factor 10 for visualization purposes. 



for z ^ 4 and by a factor 2 at z 2 — 3. In BOl the value fesc = 10% was preferred, based on comparisons of the predicted 
ionizing flux with estimates from the proximity effect at 2 ^ z ^ 4. Estimates of the UV background from th e pro ximity effect 
are known to be larger than those obtained via theoretical models of the IGM opacity (i.e. the work of Fan e t alj we are using 
here), probably because of a bias of the QSO distri bution towards th e denser environments the (iSchirber fc Bullockll20o3) 
or because of systematic errors due to line blending jScott et alJ l2000 |). Furthermo re, proximity-based UV backgrounds have 
been traditionally derived in a flat Einstein-De Sitter cosmology (Scott ct al. 2000*): this may cause a further overestimation 
by 40% if the true cosmology is the one we are using here (jPhi lliPDS ct al. 2002). 

Our conclusions on fesc are not affected by the high values of the measured photo-ionization rates for 5 ^ z ^5 6. As 
already pointed out in i\2.1\ (see Fig. there) that bump is due to a deviation of the IGM optical depth from the assumed 
smooth power-law evolution. Had we modeled in more detail the behavior of the IGM optical depth we would have also 
recovered higher values for F in the same redshift range. Therefore our statement above is valid throughout the whole redshift 
range. 

A third interesting point concerns the relative contribution of stellar and thermal emission to the metagalactic UV flux at 
high redshifts where the QSO emission drops rapidly. Although one can in principle reproduce the observed photo-ionization 
rates (thick solid line) with an escape fraction fesc ~ 2%, it is also possible that at least half of those rates are due to thermal 
emission. This is the case when a fraction rj ~ 0.5 of the feedback energy is re-radiated through thermal processes (thin solid 
line). However, thermal emission is characterized by a much harder spectrum than the stellar one. This might serve not only 
as a way to discriminate observationally between the two components, but it might have implications for the temperature 
evolution of the IGM, studied in the next Section. 

Finally, the right panel of Fig. |S| shows the He ii ionization rates. According to the plot, above 4 Ry stellar emission 
is thoroughly negligible (independent of /esc) whereas thermal emission is comparable to QSOs at 2 ~ 3 but completely 
dominates the radiation flux at higher redshifts. This result is very important in terms of the IGM evolution and has not 
previously been noticed. It depends only weakly on feedback, but it does assume an escape fraction of the thermal photons 
from collapsed halos of order 1. As already discussed in Sec. 12. 31 and shown in Appendix^ these conditions should be ensured 
for coUisionally ionized gas within halos of virial temperature above 10*' K, that is for the halos that generate most of the 
thermal emission. However, at 4 Ry thermal flux would dominate over the QSO component even for an escape fraction from 
halos ^ 0.3 at 2 ^ 4 and as small as ^ 0.05 at 2; ^ 5.5. 



3.4 He II Reionization 

The previous results hint at the intriguing possibility that He 11 reionization could have been powered by UV light from cosmic 
structure formation. For this to be the case, the production rate of ionizing photons has to satisfy the condition 
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r X mm{rrec,THubble) > 1- (14) 

Here Tree — 0.9/ngas<^{T)C' is the recombination time, a is the radiative recombination coefficient, and the clumping factor 
C = (np)/(np)^ > 1 is meant to allow for the effects of density inhomogeneities inside the ionized region. Using a helium to 
hydrogen number ratio y = 0.08 and assuming a temperature of the reionized gas T ~ 4 x 10* K, we find 

Tree = 5.3 X 10^^ C"' ("1^) ^ ^' (^^^ 

which is shorter than a Hubble time for z ^ 4.5. Thus, from Fig. |S|we find that at z ~ 6 thermal emission dominates the 
photoionization rates and alone provides 5.5 x 10~^^ (10~^^) He ii photoionizations/s in the no-feedback (feedback) case. 
According to eq. I15I I. the He ii recombination rate at the same redshift is 6.4 x 10~^'^C s~^, i.e. a comparable value. Hence, 
it appears that structure formation can produce He ii reionization around z — 6, without the contribution from any other 
process and essentially independently of the feedback prescription adopted. 

An additional check of the above results can be performed. At z — 6 the fracti on of mass in halos with virial temperature 
[calculated from Press-Schechter, eq. 0] larger than 4 Ry/fcs is equal to fi = 3% iMo fc WhitdEooS) . In order to provide at 
least one He ii ionizing photon per He atom (we have assumed a helium abundance y = 0.08 and that He is all singly ionized) , 
the required mean photon energy has to be equal to 4 Ryx (j/Z/i) = 145 eV. By averaging again over the mass distribution we 
find that this mean energy is 1.7 keV. Hence this simple argument confirms that He ii reionization can be caused by structure 
formation. 

Whether or not this possibility is fully compatible with other observational results is not clear. Nevertheless, in 331 we 
show that this scenario predicts a temperature evolution of the IGM that maybe consistent with existing observations. 



4 DISCUSSION 

Fig. |7| illustrates the thermal history of the IGM. Data points are from lSchave et all i2000h . who measured the cut-off in the 
b — Nm relation in a sample of nine quasar spectra in the redshift range 2.5-4.0. Despite the large uncertainties, there is an 
indication of a temperature peak at redshift z « 3. On the other hand, the plotted curves represents different evolutionary 
scenarios and are (numerical) solutions to the equation 



dT , 
— = (7- 



(r - A)dt _^ duH ^ _ ^ 1 



ytUukBT uh yt 7—1 



+ 7^ (16) 



where yt = ^^Vi, yt ~ ni/nn is the concentration of the i-th species, ^ = Vi^^ '^Vif^i is the mean molecular weight and 
the IGM mean density evolves as nH{z) = nH(0)(l + z)^ . The calculation accounts for photoheating from a time dependent 
ionizing background, Compton cooling on Cosmic Microwave Background photons and non-equilibrium evolution of hydrogen 
and helium ionic species in a cosmological context (see Fcrrara & GialloriEO 1996; Thcuns et al. 1998b, for further details; 
we do not include metal cooling which is negligible for typical IGM metallicities). The temperature of the medium sd, z — 10, 
where we start our integration, has been computed self-consistently. 

We emphasize that the various curves presented in Fig. |7| are meant to illustrate a variety of possibilities and not to 
depict a specific scenario. In this perspective, we first discuss the dotted line, representing the case for a UVB generated by 
QSO-|-stellar contributions given by the photoionization rates presented in Fig. 1131 In this case H i reionizes at z ~ 6 and 
He II around z ~ 3 (therefore their heating and ionization rates are set to zero prior to their respective epoch of ionization). 
At hydrogen reionization the IGM temperature jumps to log T = 4.13 and decreases later on due to adiabatic cooling. A sharp 
rise is seen at He ii reionization at ZHeii = 3, although not quite enough to fit all the data points. When the contribution from 
thermal emission is added, the effect is small since the previous photoionization rates were already sufficient to keep the gas 
fully ionized and, therefore, their increase changes little the energy input to the IGM. 

As a second example the long-dash-dot curve i llustrates a simila r case but with the photoheating rates increased twofold 
to mimic radiative transfer effects as discussed in lAbel fc Haehneld lll999l) . Again, H i reionizes at z ~ 6 and He ii around 
2 3. This scenario provides a somewhat better fit. Obviously it is possible to improve even further the agreement with the 
data points by fiddling with the a vailable parameters so as to produce the desired behavior. At least in principle the data 
points can be reproduced (see also lRicotti et aDl200(J) . 

Motivated by our findings in 13.41 we also explore a less conservative scenario in which He ii reionization occurs as early 
as zhgII = 6 and is mainly driven by thermal emission. By construction, the temperature peak around z ~ 3 does not develop 
in this case. Nevertheless, most of the data are fit at the 1 a level and all except one point at 2a level. 

The possibility that cosmic structure formation is responsible for He ii reionization at an early epoch is intriguing and 
challenges the common wisdom that such process is due to QSOs and occurs around z = 3. It is worthwhile then to revisit the 
arguments in support of the standard conjecture. These are essentially three: (i) an apparent boost of the IGM temperature 
around z = 3 which could be due to heating associated with He ii reionization; (ii) an abrupt change of the C iv /Si iv ratio 
at 2 fa 3 indicating a significant change in the ionizing spectrum; (iii) the detection of a patchy He ii Lya absorption at 
about the same redshift, suggestive of the final (overlapping) stages of reionization process. Taken together these facts seem 
to justify the conjecture above. 

However, as already pointed out bv lSchave et al.l (|200(J) . these three observations probe different structures: metal line 
ratios probe high density regions; the effective He ii Lya absorption is mainly produced by neutral gas in the voids; the IGM 
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Figure 7. IGM temperature evolution: For H I reionization redshift zhi = 6 and He II reionization redsliift ^hoII = 3 we sliow the 
optically thin case (dot) and a case in which radiative transfer effects enhance the photoheating rates by a factor 2 (dash). Finally the 
case of an early He II reionization at ^hcII = 6 is illustrated (solid). An escape fraction for stellar photons fesc = 1% is assumed. The 
data points are from Schaye et al. (2000). 



temperature measurements concern the ionized gas with density around the cosmic mean. Radiative transfer calculations 
show that the ionization fronts pro pagate first through the low den sity regions and only subsequently penetrate the denser 
clumps such as filaments and halos (lGnedij2OO0l:ICiardi et alj|2003h . 

We also emphasize the following difficulties. The existence of the peak in the IGM temper ature evolution needs further 
data to be fully confirmed. Although a wavelet analysis seems to confirm it l|Theuns et al.ll2002h . the data are also marginally 
consistent with no feature at z = 3, in which case early reionization would provide a better fit t o th e tempe rature data 
(Fig. |7|l. The con clusion about the jump in the C iv /Si iv ratio at z = 3 is still debated (e.g. iBok scnbcrg e^^^^Od; 
lBoksejiberelll99ij) . More recent studies of m etal absorption systems in QSOs spectra both with VLT/UVES (.Kim et alJl2002r 
and Keck/HIRES llBoksenberg et alJ |2003il) find no discontinuity in the C iv /Si iv ratio arou nd z = 3. On the o ther hand, we 
have shown (Sec. 3.2) that our smoothly increasing Sl is quantitatively consistent with the lHeap et al.l i200Ci) observations. 
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Finally. ISongailal il998|) points out that a considerable number of hard photons are required at 2; > 3 to reproduce the observed 
abundances of O vi . At that epoch a sizable fraction of cosmic volume has to be transparent at photon energies above the 
He 11 ionizing threshold. 

In conclusion, there is no conclusive argument against an early (2 ~ 6) start of He 11 reionization which may result from 
the hard UV field emitted during the virialization of large galaxies/small groups. Clearly, this conclusion needs to be confirmed 
by numerical simulations, although this task might be far from easy, as the correct treatment of the strongly radiating cooling 
inhomogeneities in the gas poses challenging numerical problems. Resorting to improved observations of the proximity efi^ect 
might then provide a check of the scenario proposed here. 



5 SUMMARY AND CONCLUSIONS 

We have shown that UVB ionizing photons can be copiously produced by thermal emission from shock-heated gas in collapsing 
cosmic structures. Our calculations are based on an implementation of the extended Press-Schechter theory. However, the 
estimated amount of thermal radiat ion is consistent with that inferred from an independent analysis based on observed, high 
redshift star formation rates (BOl; iLanzetta et al.ll2002^■ and the distributi on of stellar mass as a function of halo virial 
temperature as reconstructed from recent SDSS data jKauffmann et al.ll20o3) . 

Thermal radiation is characterized by a hard spectrum extending up to photon energies of order hpV ~ ksT. This is well 
above the H i and He 11 ionization thresholds for virial temperatures above 10^ K. The bulk of the emission is produced by 
halos with temperatures between 10^ K and a few x 10^ K, corresponding to masses 10^^"^=* Mq. We assume that most of 
the thermal radiation is able to freely escape into intergalactic space, which is justified for a gas that is coUisionally ionized 
and at the temperature of these halos (see Appendix El. 

We use simplified radiative transfer to compute the transmitted fiux due to QSO, stellar and thermal emissions. Impor- 
tantly, the resulting associated photoioniza tion rates, when compared to mea surements of the Lyman series Gunn-Peterson 
effect in the spectra of high redshift QSOs llFan et al.ll2002HBecker et alJl200j) . imply an escape fraction of UV ionizing pho- 
tons from galaxies, f^ac, be low a few %. This result is in agreement with very recent and independent determinations of fesc 
carried out by Fcrriandcz-S oto et al.l (2003), who set a 3a (statistical) upper limit fesc ~ 4% for a sample of spectroscopically 
identified galaxies of redshift 1.9 < 2 < 3.5 in the Hubble Deep Field. 

It turns out that near the H i ionization threshold, thermal emission is comparable to the stellar component and amounts 
to about 5-10 %, 15-30 % and 20-50 % of the total at redshifts of 3, 4.5 and higher respectively. The quoted range depends on 
the fraction of feedback energy allowed to be re-radiated through thermal processes. Near the ionization threshold for He 11 , 
the thermal contribution is much stronger. It is comparable to the QSO input already at 2 ~ 3, and it dominates for 2 > 4. 
Thus, this contribution, with a typical softness parameter 5*^ = 10 — 100, is expected to play a major role in He 11 reionization. 
In principle structure formation alone provides enough photons to produce and sustain He 11 reionization at 2 6. These 
conclusions are independent of our feedback prescriptions, which only affect low virial temperature systems. In our scenario. 
He II ionizing photons are produced primarily by relatively large collapsing structures (with T^ir ^ 10® K). 

The thermal spectrum cx u"' is very hard, with a slope in the range a « 1 — 2 depending on the importance of the 
supernova feedback we include. The latter process primarily affects the smallest systems, increasing their emissivity through 
reradiation of SN input. As these smaller objects dominate the 1-4 Ry band, the spectrum in this energy range becomes 
flatter as feedback is increased. Measuring the evolution of the thermal component of the UVB provides a powerful method 
to evaluate the importance of SNe input into the intergalactic medium. 
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APPENDIX A: OPTICAL DEPTH FOR COLLISIONALLY IONIZED GAS IN A VIRIALIZED HALO. 

Under coUisional equilibrium the fraction of neutral hydrogen and singly ionized helium are given by the ratio of t he coefficients 
for rad iative recombination and coUisional ionization. With the values for these coefficients as summarized in lTheuns et alJ 
Jl998d) we find 

^ . 1.3x10- '-^^'^ 



1.2 1 + 7^1 ^ 1-57809 ^ 
^6 1 , ^0 7 exp I — Al) 

1 + / 6.31515 \ 

~ 7 X 10 Tfi ■ ^ _^ y'o^ exp y-^^ — j (A2) 



J^HcIII J- + 6 

where we have used the notation r„ = T/IO" K. Within a virial radius 
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«.,.,.0.»„pc.,..-(i±.)-"(^)-*"(||)- 

the gas column density for absorption at 1 Ry (4 Ry) is A'' = fnfsrigaaRv, where /„ is given in ea. lAlHA2t when considering that 
H (He II ) is nearly completely ionised, /s is the hydrogen (helium) abundance by number and Ugaa ~ {^lb/^^m){pcrAc/'mp){l + 



z)^. Thus for the optical depth, r = aN, we find 



.(1 Ry) ^ fO-3 T,- (i^) exp {^) (A4) 

ri.Ry) ^ 2.5xfO-^-i±|^;(i±I)^^xp(^). (A5) 



APPENDIX B: PRESS-SCHECHTER ESTIMATE OF \{M,Z). 
The lifetime of a halo is 

rh{z, Zf, M) = THubble{z) — rHubble[Zf{M)], (Bl) 

with the formation r edshi ft, Zf{M), of a halo of mass M defined as the formation time of a progenitor with mass M/2 
jLacev fc ColdllQOl ^9^. Its relation to the running redshift z (assumed ^ 1) is statistical and M-dependent, and is 
described by 

where ujf takes values (effectively b etween a few xlO~^ and a few) according to a probability function defin ed wit hin the 
extended Press-Schechter formalism (iLacev fc Coldll99A see also approximation formulae in Kitavama & SuTol ll996l) and Sc 
and a{M) have been defined in >12.3.1I Thus, at a given redshift z the formation rate of cold gas within a halo of mass M that 
formed at Zf, averaged over the halo lifetime Th{z, Zf , M), is 

{Mcold)ume{z,Zf,M) ^ Z^""'" „ ^, + ( 1 + fi) (M, ) t.^e (B3) 
Th{z,Zf,M) 

where R (~ 1) is the mass return fraction by stars and the subscript explicitly indi cates time, as opposed to cosmic- volume. 
After introducing the characteristic star formation timescale: = Maoid/M, (e.g. IZasov|[l995l) . we rewrite eq. IIB3t as 

{M,old)u^e = (l + + ^) {M,)u,ne (B4) 

It follows that at early epochs when th <^ radiatively cooled gas is produced at a much faster ra t e than stars. However, the 
two quantities converge at late times, when rj, ^ r*. This is consistent with, e.g. . Ivan den BoschI i2002h semianalytic model 
of a Milky Way-like galaxy, in which star formation occurs at a rate proportional to the amount of available cold gas. Based 
on the diagram of Fig. 2 of his paper we readily infer that A ~ 5.5, at z ~ 9 shortly after the galaxy starts evolving, A ~ 2.0 

by 2; ~ 3 and A ~ 1.2 by 2 ~ 0. 

\ :■ of ~ Gyr or so have been inferred from observations of nearby galaxies by, e.g.. IZasovl lll995ll . although the 
author suggests higher values during the early stages of the evolution of galactic disks. Alternatively, r^, can be take n from 
semianalytic models which successfully reproduce a number of observed galactic properties. In the model offKauffman n et alJ 
il993h . the star formation rate is prescribed according to M* = OM^oid/Tdyn, where 6 — 0.1 ~ 0.3 is a control parameter and 
/vgai = 0.08^/0^/0.3 THubbie(z). lu tMs formulatiou 

r. = ^ = 0.08 f ^ ^''"""'^^^^ . (B5) 
u \ 0.3 / u 

It implies lower values than those found bv IZasovl (^9^ and, therefore, is more conservative for our present calculations. 
After using this expression in eq. I|B4|I and volume averaging over halos of a given mass M, a comparison to eq. Q leads to 



A(M, z) = 1 + J? + i /l^fliM^ - 1 + + 0-1%^ 

9 \ Th{z,zf)/ e 



l + z/ 

l + z 



3/2 



(B6) 



where the last equality holds for z 2> 1. The averaging operation is effecti vely over the halo formation time which we achieve 
through the relation in eq. IIB2II and the probability functions for ujf as in lKitavama fc Sutg (,.199&'l . In fact, for 9 ~ 0.1 — 0.3, 
values of A computed through eq. IIB6I 1 are of order of several. This is consistent with the idea that the bulk of the cold gas 
in galactic systems formed at z > 2 — 3 and, on average, h ad not yet been depl eted by star formation occurring at roughly 
constant rate before that time (see plot 2 in BOl; see also iLanzetta et alJl20ojV For the purpose of comparison in t|2.4l we 
have used 9 — 0.2 and R — 1, although results are not too sensitive to R. 
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